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SYSTEM, METHOD AND STORAGE MEDIUM FOR PREDICTING IMPACT 
PERFORMANCE OF PAINTED THERMOPLASTIC 

CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims the benefit of U.S. provisional patent 
applications serial niimber 60/234,428 filed September 21, 2000, serial number 
60/234,427 filed September 21, 2000 and serial number 60/273,648 filed March 5, 
2001, the entire contents of which are incorporated herein by reference. 

COPYRIGHT NOTICE 

[0002] A portion of the disclosure of this patent document contains material 
that is subject to copyright protection. The copyright owner has no objection to the 
facsimile reproduction by anyone of the patent document or the patent disclosure, as it 
appears in the Patent and Trademark Office patent file or records, but otherwise 
reserves all copyright protection whatsoever. 

BACKGROUND 

[0003] The disclosure relates generally to thermoplastic performance, and 
more specifically, to a method, system and storage medium for predicting impact 
performance of painted thermoplastic. 

[0004] Accurately predicting the impact performance of thermoplastic parts is 
a challenge for engineers and designers. In order to correctly predict the total 
response of the material and part to an impact event, the engineer or designer must be 
able to predict the load-displacement response of the part prior to failure and the 
failure behavior. Paint systems can significantly affect the impact performance of 
thermoplastic parts. Such paint systems can vary fi-om rigid to more flexible. In order 
to accurately predict the load-displacement response of the material prior to failure, 
the engineer must know the elastic behavior, yielding behavior, and post yield 
behavior of the painted thermoplastic material. Predicting if failure will occur, along 
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with the failure mode (e.g., ductile or brittle) and the load or displacement at failure, is 
more difficult. Nevertheless, this is necessary to determine whether a painted part will 
meet its impact specifications (typically described as energy absorption criteria). If 
the possible failure modes are not known, and if accurate failure criteria do not exist 
for each failure mode at the appropriate strain rate and temperature of the application, 
then the engineer caimot predict the energy absorption capability of the painted part. 
However, the current technique of determining impact performance by first 
manufacturing a part, and then testing the part, is wastefiil, time-consuming and 
costly. 

[0005] Finite element analysis ("FEA") is useful for predicting the structural 
performance of plastic components. Through the use of finite element tools, 
conceptual designs may be assessed and mature designs may be optimized; thereby, 
shortening the costly build and test cycle. In the past, predictions were most useful in 
predicting the load-displacement response of the component. This was done by 
accurately modeling the geometry and bomdary conditions, and by knowing the 
modulus of the material. However, as plastics are increasingly used in more 
demanding applications, such as load bearing automotive components, other nonlinear 
deformation processes and failure mechanisms become important. In plastics, the 
yield stress is typically strain rate sensitive and can be pressure dependent as well. 
Another consideration is the actual failure event (e.g., whether the material will 
behave ductiley or brittlely and under what condition will it behave ductiley or 
brittlely). 

[0006] The value most typically used to predict failure, a true strain to failure 
number, is often not measured correctly. Often a percent elongation number is 
actually reported which is not a material property, but rather depends on the geometry 
of the tensile specimen. It is simply a measure of the total elongation of the specimen 
divided by its initial length. Attempts to measure a true strain number in a tensile test 
may be difficult, since most polymers neck and locally deform. The strain needs to be 
measured locally at the point of necking, which is unknown a priori in a standard 
ASTM or ISO tensile bar. In addition, even if the point of necking is known, standard 
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extensometers are not refined enough to measure the local strain that occurs prior to 
failure. When a true strain to failure value is accurately determined using tensile data, 
it is time consuming and costly, and is usually done optically. 

[0007] Nevertheless, the failiure criterion currently used to predict whether or 
not failure will occur in an unfilled thermoplastic is typically a strain to failure value. 
Often, a percent elongation result or the equivalent is input into a finite element code. 
This is not correct, however, because a percent elongation does not represent a strain 
to failure value. The percent elongation is simply the ratio of the crosshead 
displacement at failure divided by the initial gauge length of the specimen. Since 
thermoplastic materials neck, the actual region of the specimen that is undergoing 
large deformations is smaller than the initial gage length. The displacement and 
accompanying strain is localized in the necked region. 

[0008] In addition, finite element codes require a true strain value whereas a 
percent elongation is defined as an engineering strain value. To accurately obtain a 
true strain failure value in a tensile specimen is difficult, because of the necking and 
strain localization that occurs prior to failure. Optical extensometers are used at high 
displacement rates to overcome the difficulty of mounting and holding a mechanical 
extensometer in place at high strain rates. An optical extensometer by itself is not 
sufficient, however, because the strain recorded is still measured over a prescribed 
distance and not locally at the failure point. If the specimen is gridded and the 
deformation pattern recorded optically, reasonable values of true strain can be 
obtained at the failure point. Of course these tests are more time consmning and 
costly than traditional tensile tests. 

[0009] Of even greater significance is the effect of the painted system on the 
underlying thermoplastic substrate. The painted surface layer is much more brittle 
than the underlying plastic substrate. Often brittle cracks will initiate in the painted 
material and then propagate into the underlying thermoplastic material resulting in 
less energy absorption capability. This will cause the part to fail much sooner than an 
identical unpainted part. It is important to characterize the failure performance of the 
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paint-thermoplastic system. This is sometimes done using a painted uniaxial, tensile 
specimen; however, most painted applications will see a biaxial stress state. In 
addition failure strain values obtained from painted tensile specimens often exhibit 
considerable of scatter and have not been demonstrated to be a good predictor of 
actual painted part performance. 

SUMMARY 

[0010] The above described drawbacks and deficiencies of the prior art are 
overcome or alleviated by a method for predicting impact performance of an ^icle 
constructed of a painted material. The method includes: applying physical properties 
of the material to a constitutive model; performing biaxial property tests on painted 
samples of the material shaped according to test geometries; performing finite element 
simulation analysis on the test geometries using the constitutive model; determining 
maximum principal stress levels from the finite element simulation analysis 
corresponding to experimental failure displacements obtained from the biaxial 
property tests; and applying the maximum principal stress levels and the constitutive 
model to finite element simulation analysis of the article. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0011] Referring now to the drawings wherein like elements are numbered 
alike in several FIGURES: 

[0012] FIG. 1 is a block diagram of an exemplary system for predicting impact 
performance of painted thermoplastics; 

[0013] FIG. 2 is a flow chart generally depicting a method for predicting 
impact performance of painted thennoplasics; 

[0014] FIG. 3 depicts an exemplary plot of a imiaxial stress-sfretch cxnve for 
polycarbonate; 
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[0015] FIG. 4 depicts an exemplary plot comparing analytical load 
displacement responses for a barrier impact of an automotive bumper using the von 
Mises yield criterion versus a pressure-dependent yield criterion; 

[0016] FIG. 5 depicts an exemplary plot showing yield stress of polycarbonate 
as a function of strain rate and temperature; 

[0017] FIG. 6 depicts an exemplary plot of an elastic-perfectly plastic versus 
multilinear plasticity model; 

[0018] FIG. 7 depicts an exemplary biaxial test geometry for determining 
failure criteria; 

[0019] FIG. 8 depicts an exemplary biaxial test with the paint side up; 

[0020] FIG. 9 depicts an exemplary biaxial test with the paint side down; 

[0021] FIG. 10 depicts an exemplary plot comparing analytical and 
experimental load-displacement traces from a disk; 

[0022] FIG. 1 1 depicts an exemplary plot of an equivalent plastic failure strain 
versus strain rate for a painted thermoplastic material at 70° F; 

[0023] FIG. 12 depicts an exemplary plot of maximum principal stress versus 
strain rate for a painted thermoplastic material at 70° F; 

[0024] FIG. 13 illustrates an exemplary method for obtaining parameters used 
in a deformation model; 

[0025] FIG. 14(A-E) illustrates an implicit finite element material subroutine 
for the exemplary constitutive model for modeling deformation behavior of painted 
thermoplastics; and 

[0026] FIG. 15(A-E) illustrates an explicit finite element material subroutine 
for the exemplary constitutive model for modeling the deformation behavior and 
failure behavior of painted thermoplastics. 
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DETAILED DESCRIPTION 

[0027] FIG. 1 is a block diagram of an exemplary system for predicting impact 
performance of painted thermoplastics in one embodiment. The system may include a 
host system 2, a network 4, one or more mechanical testing machines 18, one or more 
test fixtures with a target 20 and a data acquisition system 16 for use with the 
mechanical testing machines 18. One or more user systems 14 may be coupled to the 
host system 2 via the network 4. Each user system 14 may be implemented using a 
general-purpose computer executing a computer program for carrying out the process 
described herein. The network 4 may be any type of known network including a local 
area network (LAN), wide area network (WAN), global network (e.g., Internet), 
intranet, etc. Each user system 14 and the host system 2 may be connected to the 
network 4 in a wireless fashion and network 4 may be a wireless network, hi another 
embodiment, the network 4 may be the Internet and each user system 14 may execute 
a user interface application (e.g., web browser) to contact the host system 2 through 
the network 4. Alternatively, the user system 14 may be implemented using a device 
programmed primarily for accessing network 4 such as WebTV. 

[0028] The host system 2 may include one or more servers. In one 
embodiment, a network server 8 (often referred to as a web server) may communicate 
with the user systems 14. The network server 8 may be implemented using 
commercially available servers as are known in the art. The network server 8 handles 
sending and receiving information to and firom user systems 14 and can perform 
associated tasks. The host system 2 may also include a firewall server 10 to: (a) 
prevent unauthorized access to the host system 2; and (b) with respect to 
individuals/companies that are authorized access to the host system 2, enforce any 
limitations on the authorized access. For instance, a system administrator typically 
may have access to the entire system and have authority to update portions of the 
system. By contrast, a user contacting the host system 2 from a user system 14 would 
have access to use applications provided by applications server 12 but not alter the 
applications or data stored in database 6. The firewall server 10 may be implemented 
using conventional hardware and/or software as is known in the art. 
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[0029] The host system 2 may include an applications server 12. Applications 
server 12 may execute a plurality of software applications or modules as shown in 
FIG. 1. The applications may include a finite element module 30 and a design module 
40. The finite element module 30 may access a user-defined finite element material 
(UMAT) subroutine 32 and a vectorized (explicit) user-defined finite element material 
(VUMAT) subroutine 34, as will be described in further detail hereinafter. Each 
module and subroutine may serve as a tool that aids in predicting impact performance 
of thermoplastic as described herein. Note that each module may be implemented 
through a computer program. The computer program(s) that implement the modules 
may be stored on applications server 12 or may be stored in a location remote fi:om 
applications server 12. Alternatively, more than one apphcations server may be used 
to execute the software applications or modules. The finite element analysis software 
is commercially available but, alternatively, may be user specified computer code. 
For example, finite element module 30 may that which is commercially available fi:om 
Hibbit, Karlsson, & Sorensen, Inc. under the name ABAQUS. 

[0030] The applications server 12 may be coupled to a database 6. Database 6 
may contain a variety of information used by the software applications or modules. 
The database 6 may include data related to the development of a thermoplastic 
product, such as material specifications, material properties, constitutive model 
parameters, failure criteria, target 20 and part test results, and the like. The database 6 
may also include design data, such as material comparison plots, finite element 
analysis data, finite element results, data comparing design iterations and the like. 
The database 6 may be an electronic database directly coupled to the applications 
server 12, or the database 6 maybe in the form of separate electronic files, spreadsheet 
files, or the like. The data may also be stored on paper files and manually input into 
finite element analysis files. 

[0031] One or more user systems 14 and/or the host system 2 may be coupled 
to the data acquisition system 16. The data acquisition system 16 may be used as part 
of, or in conjunction with, the mechanical testing machines 18 to record the load- 
displacement response of the target 20 being tested. The data acquisition system 16 
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may record data electronically into a computer file, or other recording means, such as 
a strip chart recorder may be used. Time may also be recorded so as to check the 
displacement rate at which the test is performed. 

[0032] One or more mechanical testing machines 18 may be coupled to the 
data acquisition system 16 for testing one or more targets 20, such as a thermoplastic 
part or specimen. The mechanical testing machines 18 are used to perform material 
tests required to determine the deformation and failure behavior of the target 20 
material, as will be discussed hereinafter. The mechanical testing machines 1 8 may 
be servohydraulic machines, but other types of machines can be used. The mechanical 
testing machines 18 may be used to perform tests at displacement rates that simulate 
the strain rate(s) of interest of the apphcation or end use of the part. Also, the ability 
to test at various temperatures may be needed, unless the application of the target 20 is 
at room temperature only. Preferably, the testing is performed in an environmental 
chamber. Alternatively, the target 20 may be cooled or heated in a separate chamber 
to the temperature(s) of interest and then tested (preferably within about a minute) 
before a change in temperature. Test fixtures used for the target 20 may be used in 
conjunction with the mechanical testing machines 18 for performing various tests to 
determine characterization (e.g., tensile tests, compression tests, disk impact tests, 
notched beam tests and the like). As described later, one or more of these tests may 
be used to simulate three stress states that a part may experience in actual use: 
uniaxial, biaxial and triaxial. 

[0033] hi general, application server 12 is coded with a method for 
characterizing the aspects of material behavior into a set of material data and 
parameters needed to predict the impact performance of painted thermoplastic. A 
material transfer function (constitutive model) having five adjustable parameters that 
are determined via testing and data reduction techniques described herein is used in 
finite element analyses of the test parts to determine failure criteria. Note that the 
material transfer function simultaneously accounts for rate and pressure dependence, 
as well as ductile and brittle failvire modes. Once the constitutive model parameters 
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and the failure criteria are obtained, finite element analyses are performed by finite 
element module 30 to determine the performance of the part. 

[0034] As discussed below, the constitutive model developed by the method 
described herein includes a deformation model and failure criteria. The deformation 
model characterizes how a material will behave prior to failure (e.g., how it will 
deform in response to applied loads). The failure criteria help to identify whether the 
failure will be ductile or brittle. Further, the effects of strain rate and temperature 
upon the failure mode and failure criteria are determined. The effect of stress state on 
the failure behavior is also determined. As discussed later, three factors, stress state, 
strain rate and temperature, are considered when determining the failure behavior of a 
material. 

[0035] Using the deformation model, the load-deflection response of the 
material can be accurately predicted prior to failure. The failure criteria map out 
possible failure modes of the material based on stress state, temperature and strain 
rates. Different failure criteria are developed for the different failure modes (e.g., 
ductile and brittle). Knowing the possible failure modes of a material, along with 
having accurate failure criteria for the different failure modes, the impact performance 
of the part may be accurately predicted. The embodiment described herein will focus 
on the brittle failure mode and the determination of failure criteria for that mode. 
Testing has shown that the brittle failure of painted parts can be accurately predicted 
using the method described herein. Ductile failure criteria may be obtained using the 
method outlined in the U.S. Patent Application entitled SYSTEM, METHOD AND 
STORAGE MEDIUM FOR PREDICTING IMPACT PERFORMANCE OF 
PAINTED THERMOPLASTIC, attorney docket number 08EB03108, filed 
concurrently herewith on September 21, 2001. 

[0036] Referring to FIG. 2, a flowchart generally shows a method 100 
employed by host system 2 of FIG. 1 for predicting impact performance of painted 
thermoplastic. Method 100 begins at step 102, where deformation testing of the 
material at target 20 is performed by mechanical testing machines 18. As will be 
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described in fiirther detail hereinafter, deformation testing may include tension and 
compression testing of the material at the service conditions of the part to be 
manufactured jfrom the material. The deformation testing may be performed on a 
painted or unpainted sample. Typically, the thickness of the thermoplastic substrate is 
much greater (e.g., 30 times) than the thickness of the painted layer. As a result, the 
substrate material will dominate the deformation behavior of the system prior to 
failure. Therefore, the deformation behavior of the painted system can usually be 
determined using an unpainted sample. Method 100 continues to step 104 where a 
material deformation model is created using data jfrom the deformation testing 
performed at step 102. Preferably, step 104 includes determining pressure dependent 
material parameters by comparing tensile and compressive 3deld stresses at the same 
rate for varying temperatures. After the deformation model is created in step 104, 
method 100 continues to step 106 where failure testing is performed by mechanical 
testing machines 18 on painted material samples using biaxial property tests at the 
service conditions of the part to be manufactured from the painted material. 
Preferably, biaxial testing is performed on the painted and unpainted surfaces of the 
samples, as will be described in further detail hereinafter. Method 100 then continues 
to step 108 where finite element analysis (FEA) is performed on the test sample 
geometry from step 106 with the FEA model employing UMAT 32, which embodies 
the deformation model developed in step 104. In step 110, the deformation model of 
step 104 is validated by comparing the loads and displacement values obtained by 
fail\u:e testing in step 106 with the loads and displacement values obtained from the 
FEA model of the sample geometry obtained in step 108. If, from the comparison of 
step 110, the deformation model can be validated to a desired degree of accuracy, 
method 100 then continues to step 112. In step 112, the maximum principal sfress 
levels from the FEA analysis of step 108 are correlated with the experimental failure 
displacements determined in the failure testing of step 106. Method 100 then 
continues to step 114 where the brittle failure criteria of the material are obtained from 
the correlations of step 112. Preferably, step 114 includes using the correlations of 
step 112 to determine for each strain rate the maximum principal stress levels 
corresponding to the initiation of brittle failure, and then plotting the average 
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maximum principle stress levels as a function of strain rate for brittle failure. Finally, 
in step 116, the deformation model developed in step 104 and the failure criteria 
developed in step 1 14 are applied to FEA of the painted article to be manufactured to 
predict the impact performance of the material in the part. Preferably, the deformation 
model and failure criteria are embodied in VUMAT 34, which can be accessed by the 
FEA model. 

[0037] The general method 100 will now be discussed in further detail in the 
following sections I through IV. Section I provides general information for the 
characterization of the deformation model, as appHed in steps 102 and 104 of method 
100. Section n provides a preferred embodiment for characterization of the 
deformation model, as apphed in steps 102 and 104 of method 100. Section IE 
provides a preferred embodiment for characterization of the failure criteria, as apphed 
in steps 106-1 14 of method 100. Section IV provides a preferred embodiment for the 
application of the deformation model and failure criteria of the material in a finite 
element analysis to predict part failure, as apphed in step 1 16 of method 100. 

I. DEFORMATION MODEL CHARACTERIZATION: 

[0038] In determining the deformation behavior of a material, three areas may 
be considered; the elastic response, the yield response, and the post yield response. 

ELASTICITY: 

[0039] For thermoplastics, characterizing the elastic response of the material 
for use in finite element codes may be obtained through standard ASTM or ISO 9000 
test procedures on tensile bars. Usually an elastic modulus and Poisson's ratio is all 
that is needed. Poisson ratios for most thermoplastic resins range between about 0.35 
and about 0.4. The elastic modulus is typically not very rate sensitive, although at 
slow strain rates, somewhat lower modulus values may be obtained because of 
viscoelastic effects. Modulus is somewhat temperature dependent and, therefore, may 
be tested at the application temperature(s). 



11 



08EB03119 



[0040] Standard elastic-plastic stress-strain models in commercial finite 
element codes assume a linear elastic response prior to yielding. In actuality polymers 
exhibit a nonlinear elastic response prior to the plateau in their stress-strain curve. If 
the yielded region of the part is small, this nonlinear elastic response prior to yield can 
typically be ignored without having a noticeable affect on the global load-deflection 
prediction. If the yielded region is significant such as in simulating a dynatup 
puncture test, including this nonlinear elastic response may result in even more 
accurate predictions. This behavior is most often implemented in a user defined 
material subroutine. 

[0041] Hyperelastic material models are available in some commercial finite 
elements codes which allow for nonlinear elastic material behavior up to very large 
strains. However these models are intended for rubber like or elastomeric materials. 
These materials are essentially incompressible, do not yield, and can experience 
elongation of several hundred percent. 

YIELDING: 

[0042] Yielding is typically defined using a simple tensile test, hi plastic 
finite element simulations, the yield sti-ess is usually taken to be the initial peak in a 
uniaxial stiress-sti-ain curve (see FIG. 3). Even though this yield sti-ess limit is 
associated with a uniaxial stiess field, effective sttess expressions are available to 
define yielding for multiaxial stiess fields as well. Yield predictions can then be made 
by comparing an effective, multiaxial stress, usually the von Mises stress, given in 
equation (1), with the uniaxial yield stress of the material. Yield occurs when the 
effective, von Mises stress equals the uniaxial, tensile yield stress. Although the von 
Mises yield criterion originated for metals, it has been used successfully to predict 
load-deflection behavior in thermoplastic parts experiencing yielding. 
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a 



—f=-^{Gi - Qj)^ + (cr2 - ^3)^ + (cJi - as) =a^ 
V2 (1) 

where: 

<7mises IS the von Mises yield stress 

o-y is the uniaxial yield stress 

a^, CJ2 and 03 are the principal stresses 

PRESSURE EFFECTS UPON YIELDING: 

[0043] As a consequence of using the von Mises criterion, yielding is assumed 
to be independent of hydrostatic stress or pressure defined in equation (2): 



(2) 



where: 



(7h is the hydrostatic stress 

, and cfg are the principal stresses 

[0044] However, many thermoplastics do display pressure-dependent yielding 
behavior. Tensile hydrostatic stresses tend to decrease the yield stress, while 
compressive hydrostatic stresses tend to increase the yield stress. Note that the 
hydrostatic pressure is equal in magnitude to the hydrostatic stress but opposite in 
sign, (e.g., multiplying the hydrostatic stress by negative one gives the hydrostatic 
pressure). 
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[0045] In most cases, pressure effects on yielding are not significant from an 
engineering viewpoint and may be ignored. Since most thermoplastic parts are thin 
walled, large hydrostatic stress fields cannot develop, except possibly near some local 
stress concentrations. Therefore, ignoring pressure effects will not significantly affect 
gross part performance predictions. For example, a comparison of load-deflection 
behavior for a barrier impact of an automotive bumper using a von Mises yield 
criterion versus using a pressure-dependent yield criterion is shown in FIG. 4. Note 
that little difference is observed when a pressure-dependence parameter characteristic 
of a Polycarbonate/Polyester blend is used. 

[0046] For materials with a large rubber content, pressure effects on yielding 
are more significant. For these materials, cavitation of the rubber occurs under tensile 
stress fields resulting in a lower tensile yield stress. In a standard tensile test these 
materials typically experience large extensions (e.g., greater than about 50%) with 
little or no lateral contraction because of the "additional volume" created by the 
cavitation of the rubber. Under compressive stress fields the rubber does not cavitate 
resulting in a larger compressive yield stress. If a material's yield stress is 
significantly pressure dependent, and if a part sees large regions of compressive stress, 
then using a pressure dependent yielding model (or a separate tensile and compressive 
yield stress) in a finite element analysis would yield better part performance 
predictions. A comparison of the tensile and compressive yield stress of a variety of 
polymers at room temperature is shown in Table 1 . 



14 



08EB03119 



TABLE 1 



Material 


Tension (mega 
Pascals) 


Compression 
(mega Pascals) 


Ratio 


Acrylonitrile-butadiene-styrene 


42.5 


45.9 


1.08 


Polycarbonate/Acrylonitrile- 
butadiene-styrene 


56.0 


62.0 


1.11 


Polycarbonate 


66.0 


66.0 


1.00 


Noryl® GTX910 


55.0 


76.0 


1.38 


Noryl® EM6100 


36.0 


60.0 


1.67 


Polybutylene terephthalate 


53.1 


73.1 


1.38 


Polycarbonate/Polyester 


51.0 


61.0 


1.19 



STRAIN RATE AND TEMPERATURE EFFECTS UPON YIELDING: 

[0047] The yield stress of a polymer depends upon the rate and temperature at 
which it is tested. In general, higher rates and lower temperatures lead to higher yield 
stresses. Examples of temperature and rate effects on the yield stress of polycarbonate 
are shown in FIG. 5. As shown, the yield stress increases linearly for each order-of- 
magnitude increase in strain rate. 

MATERIAL CHARACTERIZATION AND MODELING FOR YIELDING: 

[0048] When yielding is included in a numerical simulation, yield data at the 
appropriate temperatures and strain rates is used. The temperature experienced by the 
part is usually known. However, the strain rate experienced by the part is calculated. 
The strain rate may be approximated by dividing the maximimi strain found in the part 
at a given displacement by the time it took the part to reach that displacement. If the 
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part geometry and loading is simple enough, the strain may be calculated using 
closed-form solutions. For more complex geometries and loadings, or for more 
accurate results, an elastic fmite-element analysis may be performed to calculate the 
strains in the part for a given deflection. 

[0049] Once the temperature and strain rates of the part are known, tensile 
testing maybe performed under these conditions to determine the temperature and rate 
dependence of the yield stress. Note that an application is usually at a imiform 
temperature when being impacted, therefore simulations may be performed without a 
temperature dependent yielding model (provided the material has been tested at the 
application temperature). For example, if an apphcation must meet certain energy 
requirements at room temperature and -30° C, then tensile testing may be performed 
at each temperature (with the appropriate data being used to simulate performance 
under each temperature). Estimating the strain rate in a component is more difficult. 
In addition, the strain rate in the part will vary from location to location. 

[0050] hi certain finite element analyses that account for dynamic effects, 
strain rates are calculated internally and a rate-dependent yielding model may be 
defined. This approach eliminates the need for performing an initial elastic finite 
element analysis to estimate strain rate and allows the yield stress to vary throughout 
the part based on local strain rates. If a rate dependent plasticity model is used, testing 
may be performed over a range of strain rates and a rate dependent plasticity model fit 
to the data. For most polymers, yield stress varies linearly versus the log of strain rate. 
Preferably, yield stress values are measured over a few decades of strain rate and 
cover the range of strain rates encountered in the apphcation. 

[0051] If a rate dependent yielding model is not going to be used, the strain 
rate experienced by the part is estimated. Since yielding will occur first in the most 
highly strained region, it is recommended that the strain rate be calculated for the 
region of highest strain. Since the effect of rate on yield stress is only significant for 
orders-of-magnitude variations m rate, approximating the yield stress using the 
maximum strain rate in the part is usually sufficient. If the strain rate is beyond 
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testing limits, the yield stress may be tested over a few orders of magnitude of strain 
rate. A linear fit of yield stress versus log strain rate can be used to extrapolate the 
yield stress out to the higher strain rate. 

POST YIELD BEHAVIOR: 

[0052] For many polymers, a decrease in stress is seen immediately after yield 
followed by a subsequent increase in stress (see FIG. 3). These post yield behaviors 
are referred to as strain softening and strain hardening, respectively. This hardening 
behavior is caused by molecular chain alignment. Post yield behavior may be 
important for predicting the performance of structures experiencing areas of high 
strain. For many thermoplastics, strain hardening occurs for strains beyond about 
40%. If strains are expected to be below about 40%, the simplest model to use in a 
finite element analysis is an elastic-perfectly plastic model. In this model, an elastic 
modulus and yield stress are entered. The stress-strain curve is assumed to remain flat 
following yield, e.g., perfectly plastic (see FIG. 6). For strain levels larger than about 
40%, most polymers begin to display strain hardening behavior. In such a case, a 
muhilinear plasticity model (accounting for strain softening, and more importantly for 
strain hardening behavior) will provide more accurate results. 

MATERIAL CHARACTERIZATION FOR POST YIELD BEHAVIOR: 

[0053] To determine the yield stress of the material, it is sufficient to perform 
a standard tensile test and measure the engineering stress-strain response (which is 
based on the initial cross sectional area and gage length of the specimen) as shown in 
equations (3) and (4). This results in a stress-strain curve that rises to reach a peak 
and gradually decrease. 

o = F/Ao (3) engineering stress 

€ = (l-lo)/lo (4) engineering strain 
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where: 

a is the engineering stress 

s is the engineering strain 

F is the load on the specimen 

Ao is the initial cross sectional area 

lo is the initial gage length 

1 is the current gage length 

[0054] If the post yield behavior of the material is desired, then a true stress- 
strain curve is used. A true stress-strain curve is more difficult to obtain since the 
stress is based on the current cross sectional area, and not the initial cross sectional 
area. Once necking initiates, the cross sectional area changes quickly, resulting in an 
initial load drop. Then the material starts to harden and the neck propagates, resulting 
in a increase in the true stress response, (e.g., strain hardening). True stress and strain 
equations are shown in equations (5) and (6). 

fft = F/Ai (5) true stress 

£t = ln(l/lo) (6) true strain 

where: 

at is the true stress 

St is the true strain 

F is the load on the specimen 

Ai is the instantaneous cross sectional area 

lo is the initial gage length 
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1 is the current gage length 

[0055] Two general measurement techniques may be used for measuring the 
true stress-strain response. One option is to grid the specimen and optically record the 
deformation while recording the load, as is known in the prior art. The deformation 
measurements may be made in the necked region. This technique is accurate, 
although tedious. Another technique is to perform compression testing on cylindrical 
specimens. In a compression test the measurement difficulties associated with 
necking are eliminated since the specimen diameter expands uniformly. Once 
plasticity occurs, the material behaves incompressibly, e.g., it is volume preserving. If 
volimie is preserved, relationships relating true stress and strain to engineering stress 
and strain are obtained (see equations (7) and (8)). Compressive engineering stress- 
strain data may be recorded and converted to true-stress stain data through equations 
(7) and (8). 

at= a(l+s) (7) true stress to engineering stress 
St=bi(l+s) (8) true strain to engineering strain 

where: 

at is the true stress 

St is the true strain 

a is the engineering stress 

s is the engineering strain 

[0056] When performing a compression test, some practical considerations are 
noted. First, certain specimen dimensions may be important: specifically the ratio of 
the height of the specimen to the diameter of the specimen. If the height to diameter 
ratio is too large, buckling may occur. If the height to diameter ratio is too small 
barreling may occur. Barreling refers to the specimen taking a "barrel" shape, e.g., its 
sides bulge out at the center. Barreling is caused by firictional forces restraining the 
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lateral growth of the specimen where it contacts the platen. This may resuh in a non- 
uniform cross section as well as "dead" conical zones adjacent to the platens where no 
deformation occurs causing erroneous stress-strain data. Barreling may be minimized 
by reducing fidction between the specimen and the platen through lubricant sprays or 
sheets. To minimize barreling without buckling, a length to diameter ratio of about 
1.5 to about 2.0 is recommended, hi addition, compression tests may be performed at 
low strain rates to avoid heating the specimen (which would soften the stress-strain 
response). Note that strain rates on the order of 1x10"^ to 1x10'^ 1/s may be tested 
without specimen heating. 

[0057] Although using compression data to estimate the post yield behavior of 
the material may need to include practical considerations, the test may be much 
simpler to perform than the optical tensile technique. Therefore, if the rate 
dependence of the hardening is deemed important, then the optical technique may be 
preferred. On the other hand, if speed and costs are more important, then the 
compression test may be preferred. 

[0058] If the pressure dependence of the yield stress is desired, then tension 
and compression tests may be performed. The pressure dependent yielding parameter 
can be calculated by using tensile and compressive yield stress values at the same 
strain rate. The exact calculation of the pressure dependent yielding parameter will 
depend on the pressure dependent model employed. 

n. PREFERRED METHOD OF DEFORMATION MODEL 
CHARACTERIZATION: 

[0059] This section provides a preferred method of deformation model 
characterization, as apphed in steps 102 and 104 of method 100. The following 
embodiment includes a preferred technique, along with alternative techniques. This 
constitutive model may be in the form of computer code and implemented as a user 
defined material subroutine for use with standard finite element codes, as shown in 
FIG. 1 as UMAT 32 and as provided in FIG. 14. 
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[0060] To characterize the deformation behavior of a material, tension and 
compression tests may be performed on unpainted samples of the material at the 
temperatures of interest (steps 200 and 202 of FIG. 13). Tensile tests (step 200) may 
be performed on standard ASTM or ISO bars at three to four displacement rates 
covering three to four orders of magnitude in displacement rate. Typically, specimens 
may be tested at 3.81, 38.1 and 381 millimeters per second. Tests may be performed 
at higher strain rates if desired. Displacement rates are chosen to cover the range of 
strain rates that may be seen in the application utihzing the material. Compression 
tests may be performed on cylindrical specimens, with a height to diameter ratio 
between 1 to 1 and 2 to 1. Typically, specimens of about 6.2 millimeters high by 
about 6.3 millimeters in diameter may be tested. Note that if the height to diameter 
ratio is too large (e.g., greater than about 2 to 1), buckling of the specimen may occur. 
On the other hand, if the height to diameter ratio is too small (e.g., less than about 1 to 
1), barreling may occur. Note that barreling refers to the specimen taking a "barrel" 
shape. Barreling is caused by frictional forces restraining the lateral growth of the 
specimen where it contacts the platen. This results in a nonuniform cross section as 
well as "dead" conical zones adjacent to the platens where no deformation occurs, 
thus, causing erroneous stress-strain data. Teflon® coating may be applied between 
the top and bottom surfaces of the compression specimens and the platen surfaces to 
reduce friction between the two surfaces and minimize barreling. Compression 
specimens may be tested at slow rates to avoid heating of the specimens during 
testing. To avoid heating, displacement rates on the order of 0.007 millimeters per 
second may be used for the compression tests. In addition, for low temperature tests, 
tensile and compression specimens may be tested in an environmental chamber and 
allowed to equilibrate in the chamber for at least one hour prior to testing. 

[0061] Tensile tests may be performed for two primary purposes, to determine 
the elastic modulus of the material and to determine the strain rate dependence of the 
yield stress. Tensile tests may be performed in accordance to ASTM D638. 
Appropriate ISO 9000 standards may be substituted. Elastic moduli may be recorded 
in accordance with these standards. A Poisson ratio of 0.4 may be assumed for most 
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thermoplastics, with typical values ranging between 0.35 and 0.42. Note that the load 
deflection response of the material is not sensitive to this parameter, but it may be 
tested for if desired. For each displacement rate tested the yield stress may be 
recorded, with the yield stress taken as the initial peak in the engineering stress-strain 
curve. Next, yield stress may be plotted versus the natural log of strain rate. 
Typically, the yield stress varies linearly versus the log of strain rate for most 
thermoplastics. A linear regression of the natviral log of strain rate versus yield stress 
maybe obtained to characterize the rate dependence of the yield stress. 

[0062] After yielding, most thermoplastics display strain-hardening behavior, 
which may be preceded by some initial strain-softening behavior. This behavior is not 
characterized by a traditional tensile engineering stress- strain curve, which will 
flatten out and drop after yielding (since the change in the cross sectional area at the 
necked region is not accounted for). In order to characterize the post yield behavior, a 
true stress-strain curve is used. A true stress-strain cxirve is more difficult to obtain in 
tension since the stress is based on the current cross sectional area and not the initial 
cross sectional area. Once necking initiates, the cross sectional area changes quickly 
resulting in an initial load drop. Then the material starts to harden and the neck 
propagates, resulting in an increase in the true stress response, e.g., strain hardening. 
Thus, a compression test may be preferred as an altemative to performing tensile tests 
to obtaining post yield true stress-strain behavior. From a practical viewpoint, the 
compression test is a simpler test to perform, requu-es less specialized equipment 
(such as optical measurement devices), and is quicker and less costly. 

[0063] A compression test at each temperature of interest may be performed to 
characterize the true stress-strain, post yield, behavior of the material. Preferably, 
about 5 specimens at each temperatures of interest may be tested (to account for 
variation). Since necking does not occur in compression, a near uniform expansion in 
cross sectional diameter can be obtained if barreling is minimized. Note that barreling 
can be minimized by keeping a specimen height to diameter ratio about 2 to 1 and by 
reducing friction between the top and bottom surfaces of the specimen and the test 
platen. The change in diameter can be accurately predicted by assuming 
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incompressible material behavior after necking, e.g., a Poisson's ratio of 0.5. 
Assuming incompressibility, e.g., the material is volume preserving, true stress and 
true strain values can be calculated by knowing the corresponding engineering stress 
strain values through equations (7) and (8) as previously described. 

[0064] When using compression data to estimate the post yield behavior, there 
are factors that need to be accounted for. First, the compressive and tensile yield 
stresses may be different at the same strain rate. Second, the compressive test may 
need to be performed at a much slower strain rate than the part will experience. Given 
these circumstances, one technique is to use the tensile yield stress at the appropriate 
strain rate and superimpose the post yield behavior. This may be done by making a 
table of Act versus strain by picking off stress-strain points from the compression 
stress-strain curve and subtracting off the compressive yield stress. To calculate the 
total stress to be entered into the finite element code, add the Ao value from the table 
to the rate and temperature dependent tensile yield sfress, as shown in equation (9). 
Note that when using this technique, it is assumed that the post yield behavior is not 
rate dependent. 

<^totai =o-(^,r),,,,,+Aa (9) 

where: 

atotai is the true total sfress 

cr(s,T)yi^ij is the strain rate, s , and temperature, T , dependent tensile yield stress 
Act is the total stress in compression minus the compressive yield sfress 

[0065] Post yield, true sfress-sfrain curves can also be generated using tensile 
specimens, if the specimen is gridded and the deformation is recorded optically so that 
the change in cross sectional area at the neck can be accurately measured. Since these 
measurements are difficult, the necked region may be cut from the sample and the test 
continued for the propagated-necked region to determine hardening behavior. 
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However, the optical approach may be more difficult and time consuming, and 
possibly less accurate than compression tests. 

[0066] A compression test may also be performed to determine the pressure 
dependence of the yield stress at each temperature of interest. Unlike metals, most 
thermoplastics display some pressure dependent material behavior. Yielding is not 
independent of the hydrostatic stress or pressure as is assumed when a standard von 
Mises yield criterion is utilized. Tensile hydrostatic stresses tend to decrease the yield 
stress, while compressive hydrostatic stresses tend to increase the yield stress. For 
thermoplastic materials with a large rubber content, pressure effects on yielding are 
more significant. For such materials, cavitation of the rubber occurs under tensile 
stress fields resulting in a lower tensile yield stress. Under compressive stress fields, 
the rubber does not cavitate resulting in a larger compressive yield stress. By 
determining the tensile and compressive yield stresses at the same strain rate, a 
pressure dependent material parameter may be calculated and utilized in the 
constitutive model that may be represented by a user defined, finite element material 
subroutine (UMAT 32 of FIG. 1). An example of a user defined finite element 
material subroutine is provided in FIG. 14. 

THE CONSTITUTIVE MODEL: 

[0067] As previously mentioned, the constitutive model described herein 
accounts for both rate and pressure dependent plastiticity. Post yield strain softening 
and strain hardening are also accounted for. The constitutive model is shown in 
equation (10): 

6^, = Qxip[A(T){c7 - S(£^,)}] X oxp[-paA(T)] 

(10) 

where: 

is the equivalent plastic strain rate 
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Spi is the equivalent plastic strain 

A, Sq are rate dependent yield stress parameters which depend on temperature (7) 

(note that A and Sq are described in step 204 below) 

a is the equivalent von Mises stress 

S is internal resistance stress (post yield behavior) 

a is pressure dependent yield stress parameter 

[0068] A standard isotropic elasticity model is employed to model elastic 
behavior prior to yield. The elastic parameters input into the model include an Elastic 
modulus and a Poisson's ratio. 

[0069] Note that five material parameters are needed for use with the 
constitutive model as well as a post yield stress-strain table. These 5 material 
parameters are constants for a given material. Each material would have a unique set 
of constants for a given temperature. When these constants are used in the 
constitutive model, which is a mathematical representation of material stress-strain 
behavior, the stress -strain behavior of the material can be predicted which would 
enable one to predict the load-defLection response of a part, the stiffiaess of a part 
,when yielding will occur in a part, and stresses and strains in a part. The post yield 
stress-strain table describes the true stress-strain behavior of the material after yielding 
has initiated. Each material would have a unique post yield stress-strain table for a 
given temperature. 

The five parameters are shown below: 
E is the elastic modulus 
u is the Poisson's ratio 

A, are rate dependent yield stress parameters which depend on temperature (T) 
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a is the pressure dependent yield stress parameter 

[0070] Referring to FIG. 13, an embodiment of steps 102 and 104 of method 
100 is shown in further detail. FIG. 13 is an exemplary embodiment for obtaining the 
five parameters used in the constitutive model as well as the post yield true stress- 
strain table will now be described. 

[0071] In step 200, tensile tests may be performed at a temperature of interest 
over a range of rates. Tests must be performed at a minimum of two strain rates. 
Three or more strain rates are preferred for accuracy covering three or more decades 
of strain rates. Note that it may be preferable to choose strain rates that match strain 
rates typically seen in the application utilizing a particular thermoplastic part. Thus, 
about three to five replicates may be tested at each rate and temperature combination 
(to account for variation). Load displacement data may be collected and converted to 
stress strain data using standard ASTM or ISO9000 procedures. An elastic modulus, 
E, and Poisson's ratio, u, may also be calculated using standard ASTM or ISO 9000 
procedures. The yield stress is recorded for each strain rate tested, with the yield 
stress taken as the initial peak in the engineering stress-strain curve. Next, a plot of 
the natural log of strain rate vs. yield stress may be generated for each temperature of 
interest. The natural log of strain rate may be plotted on the y-axis and yield stress 
plotted on the x-axis. Typically, the yield stress varies linearly versus the log of strain 
rate for most thermoplastics. A linear regression of the natural log of strain rate 
versus yield stress is obtained and the slope, m, and y-intercept, b, of the plot is 
recorded and/or saved. 

[0072] In step 202, compression tests may be performed at a temperature of 
interest. Preferably a single slow strain rate on the order of about 0.0001 to 0.001 1/s 
may be chosen to help avoid material heat up during the test. Note that the technique 
for performing the compression test was previously described. Preferably, five 
replicates may be performed at each rate (again, to account for variation). Load 
displacement data is collected and converted to true stress-strain curves using the 
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procedure previously described. For each temperature tested, the initial peak in the 
true stress-strain curve is recorded and/or saved as the compressive tensile yield stress. 

[0073] In step 204, the pressure dependent yield stress parameter, a, is 
calculated by comparing the compressive and tensile yield stress values at the same 
strain rate. Typically, the compressive yield stress value at the low strain rate is 
extrapolated to the faster strain rate value of the tensile yield stress tests by using the 
slope, m, of the natural log strain rate versus tensile yield stress plot (as previously 
mentioned in the tensile testing description). The pressure dependent yield stress 
parameter, a, may be calculated using equation (1 1) shown below: 

3(0-; -o-;) 
"=7^>^ '''' 



where: 

a is the pressiire dependent yield stress parameter 
Cfy is the compressive yield stress 



<y^y is the tensile yield stress 

[0074] In step 206, the two parameters, and A, which determine the rate 
dependence of the yield stress, may be calculated knowing a and the slope, m, and y- 
intercept, b, of the natural log of strain rate versus tensile yield stress graph (as 
previously described). 




27 



08EB03119 



where: 

A, ^0 are rate dependent yield stress parameters which depend on temperatvire (7) 

m is the slope of the natural log strain rate versus tensile yield stress plot 

b is the y-intercept of the natural log strain rate versus tensile yield stress plot 

a is the pressure dependent yield stress parameter 

e is the base of the natural log; e = 2.71828 

[0075] Having the 5 material parameters for the constitutive model, E, u, a, 
^0 and A, the final data to be determined is the post yield behavior. 

[0076] hi step 210 the post yield stress-strain table is obtained from the 
compression data. The post yield behavior is obtained from the true stress-strain 
compression data previously obtained. A table of Aa versus plastic strain is generated 
by selecting stress-strain points from the compression stress-strain curve. Points are 
selecting starting at the initial peak in the stress strain curve, which corresponds to the 
compressive yield stress. A minimum of three points is required to define the post 
yield behavior; the yield point, the point where strain hardening begins and the end 
point of the test. More data pairs are recommended to better represent the shape of the 
stress- strain curve. Usually 8-10 points are selected are roughly evenly spaced 
intervals of strain of approximately 0.1 in/in. More or fewer points could be chosen. 
The Aa value for the table is calculated by subtracting off the compressive yield stress 
using equation (12) shown below from the true total stress value: 

Acr = a,„,,,-a;(T) (12) 

where: 

CTtotai is the true total stress 

a" (T) is the compressive yield stress at the temperature of interest 
28 



08EB03119 



Ao" is the total stress in compression minus the compressive yield stress 

[0077] The Ac versus plastic strain table is then input into the user defined 
finite element subroutine along with the five constitutive model parameters. As 
previously noted, an exemplary embodiment of the user defined finite element 
material subroutine is provided in FIG. 14. Thus, in step 208, the constitutive model 
is completely defined. 

m. FAILURE CRITERIA CHARACTERIZATION 

[0078] When assessing impact performance, failure is a concern. Most 
impacted thermoplastic parts are specified to absorb certain impact energy without 
failing. Automotive bumper impacts and impacts of electronic enclosures dropped 
fi-om height are good examples. Furthermore, two different failure modes are 
possible: ductile and brittle. 

[0079] In a ductile failure, the part fails in a slow, noncatastrophic manner in 
which additional energy is required to further spread the damage zone, hi contrast, a 
brittle failure is characterized by a sudden and complete failure that, once initiated, 
requires no further energy to propagate. Note that the failure criteria for the two 
failure modes differ. Generally, effective stress (von Mises stress) is used to assess 
when plastic (permanent) deformation has initiated. If some permanent deformation 
is acceptable, then a strain-to-failure criterion may be used as the ductile failure 
criterion indicating when tearing is expected to occur. For a brittle failure criterion, 
maximum principal stress is used to assess failure and predict part performance. 

[0080] Referring again to FIG. 2, steps 106-114 of method 100 will now be 
described in further detail. In step 106, failure testing is performed using a biaxial 
failure test. Preferably, standard Dynatup® test samples having a disk geometry, as 
shown in FIG. 7, are performed. Standard disks as used in a Dynatup® test (e.g., 
about 100 millimeters in diameter and about 3 millimeters thick), may be tested in a 
servohydraulic machine at constant displacement rates. Prior to testing, the samples 
are painted, preferably using the same paint system that will be used for the 
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manufactured part. Table 2 provides an exemplary description of two different paint 
systems that may be used. 

TABLE 2 





Primer / 


Basecoat 


Clearcoat 




Surfacer 






Paint 


Conductive, 


Water-bom 


2-component 


System 1 


solvent-bom 
flexible 
primer 


base-coat 


urethane 


Film Build 


(0.9-1.1, 
mils) 


(1.0 -1.2, mils) 


(1.8 -2.2, mils) 


Paint 


Non- 


Solvent-bom 


1 -component. 


System 2 


conductive, 
solvent-bom, 
rigid primer 


base-coat 


solvent-bom 


Film Build 


(0.8-1.0, 
mils) 


(1.0 -1.2, mils) 


(1.6 -1.8, mils) 



[0081] The following is an exemplary method for making a test sample for 
Paint System 1 : 

1 . Power wash substrate, 1 5 minute dry-off at 93° C 

2. Allow panels to cool to ambient temperature 

3 . Apply primer to achieve 0.9 - 1 . 1 mil film build 

4. Flash for 3 minutes, bake for 23 minutes at 1 1 5° C 
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5. Apply basecoat to achieve 1 .0 - 1 .2 mil film build 

6. Flash for 5 minutes at 65° C 

7. Allow panels to cool to ambient temperatures 

8. Apply clearcoat to achieve 1 .8 - 2.2 mils fihn build 

9. Flash for 1 0 minutes, bake for 23 minutes at 115 °C 

[0082] The following is an exemplary method for making a test sample for 
Paint System 2: 

1 . Wipe substrate with isopropyl alcohol and allow to dry 

2. Simulated E-coat for 40 minutes at 204 °C 

3 . Apply primer to achieve 0.8 to 1 .0 mil film build 

4. Flash for 15 min, bake for 30 minutes at 160 °C 

5. Apply basecoat to achieve 1 .0 to 1.2 mils film build 

6. Flash for 5 minutes, apply wet-on- wet topcoat 

7. Apply topcoat to achieve 1 .6 to 1 .8 mils fihn build 

8. Flash for 15 minutes, bake for 30 minutes at 140 °C 

It will be recognized that these paint systems are exemplary, and that the preferred 
paint system used will be that which will be used in the part to be manufactured. 

[0083] Disk specimens may be clamped in a rigid fixture with a clamping 
diameter of about 7.62-centimeter (3-inch) and impacted by a metal, hemispherical 
impact head with a diameter of about 12.5 millimeters. In a first series of tests, disk 
specimens having one painted surface are tested with the impact head impacting on 
the painted surface (paint side up), as shown in FIG. 8. In a second series of tests, 
disk specimens having one painted surface are tested with the impact head impacting 
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on the unpainted surface (paint side down), as shown in FIG. 9. These tests produce 
biaxial stress states. For each test, the load displacement trace is recorded and the 
displacement at break was recorded along with the failure mode, either ductile or 
brittle. Preferably, three or more specimens are tested for each condition. 

[0084] For each geometry, tests may be performed at the temperature of 
interest (application temperatures) over a range of strain rates. A small range of strain 
rates may be chosen which brackets the strain rates expected to be seen in a specific 
part or apphcation, or a large range of strain rates may be chosen, covering a few 
orders of magnitude, to more fully characterize a material. To estimate a strain rate 
that a part may experience, the impact velocity or displacement rate is determined. 
The impact velocity of interest rhay be predefined as in a regulatory, agency, or 
manufacturer required part test, or may be calculated fi-om boimdary conditions as 
specified. A common impact test is a drop test. In a drop test, a part is dropped from 
a known height or, alternatively, an object at a known height is dropped on to a part. 
The velocity at impact is calculated by equating the initial potential energy prior to the 
drop to the kinetic energy just before impact. The impact velocity may be calculated 
using equation (13) shown below: 

v = V2jA (13) 

where: 

V is the impact velocity 

g is the gravitational constant 

h is the drop height 

[0085] The strain rate may then be approximated by using closed form 
solutions, or by performing finite element analyses. Note that if a part already exists, 
the strain may be determined by instrumenting the part with strain gages, or by using 
other strain measuring techniques on the actual part. However, such techniques 
requiring an actual part are time consuming and costly. For more complex geometries 
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and loadings, or for more accurate results, an initial, elastic, finite element analysis 
may be used to estimate the application strain rate. Since the application temperatures 
are generally known, tests may be performed at those temperatures. Note that if the 
application temperature is not known, a test may be performed, preferably, at the 
coldest temperature that the part is expected to experience. Additionally, a test may 
be performed at room temperature. The test specimens may be tested in an 
environmental chamber and allowed to equilibrate at appHcation temperature for at 
least one hour prior to testing. For each test performed, load-displacement data may 
be recorded, including the displacement at break. The displacement at break may be 
recorded for each test to determine stress and strain levels at failure. The failure mode 
of the specimen, ductile or brittle, maybe recorded as well. 

[0086] hi step 108, finite element analyses are performed on the test geometry 
(e.g., the standard Dynatup® disk geometry) using the deformation model described 
earlier as input to the finite element model. The finite element model is then vahdated 
(step 110) by comparing the analytical load-displacement response from the finite 
element analyses to the experimental load-displacement results obtained fi-om failure 
testing of step 106. For example, FIG. 10 is a plot comparing the analytical load- 
displacement response to the experimental load-displacement results for an exemplary 
disk impact. Verification of the fmite element model is attained when the analytical 
load-displacement response closely resembles the experimental load-displacement 
results, within a desired degree of accuracy, hi the example of FIG. 10, for a given 
load, the experimental displacement will be generally within the range of +20% and - 
0% of the analytical displacement. 

[0087] Li the example of FIG. 10, five tests were performed. For each strain 
rate and temperature of interest, about five to ten replicate specimens are preferred (to 
account for variation). If a specimen starts to transition from a ductile-to-brittle mode, 
or if a specimen fails brittlely, then about ten to twenty tests are preferred (to account 
for the higher variability often seen in brittle failures). Also, if a significant amount of 
scatter is seen in the results then ten to twenty tests are preferred. Again, note that the 
exact number of specimens may be left to the judgment of the designer, and/or any 

33 



08EB03119 



guidelines in use and/or any statistical analysis methods that may be employed. 

[0088] After the finite element model is validated, it is then used along with 
experimental results from failure testing of step 106 to determine ductile and brittle 
failure criteria. 

[0089] Ductile failures are characterized by a tearing type event typically 
characterized by a local strain-to-failure. Since tearing failures are localized, test 
procedures to determine the true local strain at failure are difficult to define. To 
overcome this measurement difficulty, strain-to-failure values may be determined by 
correlating mechanical test results to detailed finite element analyses of the simple test 
geometries described previously. An equivalent plastic strain-to-failure is widely 
accepted as a ductile failure criterion in the art, and may be used in commercially 
available finite element packages. 

[0090] Equivalent plastic strain-to-failure values (e.g., peak equivalent plastic 
strain levels corresponding to the experimental failure displacements) may be 
obtained from the finite element predictions. For example, in FIG. 11, equivalent 
plastic strain-to-failure values are plotted as a fimction of sfrain rate for each 
temperature. The determination of ductile failure criteria is provided in detail in the 
U.S. Patent Apphcation entitled SYSTEM, METHOD AND STORAGE MEDIUM 
FOR PREDICTING IMPACT PERFORMANCE OF PAINTED THERMOPLASTIC, 
attorney docket number 08EB03108, filed concurrently herewith on September 21, 
2001. 

[0091] Brittle failures are characterized by a fast fracture usually resulting in 
specimens or parts that are broken into a few pieces, or many separate pieces. A 
brittle failure criterion, in the form of a rate dependent, critical, maximum principal 
stress criterion is used. For example, brittle failure occurs when the maximum 
principal stress in the part reaches a critical, rate-dependent, value. If maximum 
principal stress levels within the part are kept below these critical values, then brittle 
failure is not a concern. To determine critical maximum principal sfress levels that 
may initiate brittle failure, finite element analyses are performed on each test 
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geometry that failed brittley (using the deformation model described previously). 

[0092] In step 112, the maximum principle stress levels from the finite 
element analysis of step 108 are correlated with the experimental failure 
displacements determined in failure testing of step 106. hi other words, maximum 
principle stress levels corresponding to the experimental failure displacements are 
determmed using the validated finite element model. The maximum principal stress 
levels are determined for each strain rate tested, and at each temperature, 

[0093] The maximum principal stress levels for a given strain rate and 
temperature should be consistent across paint orientation (i.e., it should be consistent 
across paint side up and paint side down testing). For example, the maximum 
principal stress level predicted for each paint orientation should be the same or nearly 
the same. Preferably, values within about 20% of each other may be acceptable to 
predict part performance. Of course, a different level of consistency may be chosen. 
For example, conservatively, lower values may be chosen, and/or standard deviations 
may be calculated (probabilities of failure may be calculated given variation in failure 
criteria and/or part operating conditions). 

[0094] After the maximum principal stress levels are determined for each 
strain rate and temperature, the brittle failure criteria of the material is determined 
(Step 1 14) by taking an average of the maximum principal stress levels for each strain 
rate and temperature. The average maximum principal stress level is then plotted as a 
fimction of strain rate for each temperature. An example of such a plot is shown in 
FIG. 12. Alternatively, the lower-bound of the maximum principal stress levels for 
each strain rate and temperature may be used in lieu of the average. If a statistically 
significant number of samples have been tested, then the maximum principal stress 
may be treated statistically, thereby, establishing means and standard deviations. 
Preferably, from about ten to about twenty specimens may be tested for each set of test 
conditions. Fewer specimens, for example, about five, may be tested if the scatter in 
the failure displacement is considered low. Note that the exact number of specimens 
may be left to the judgment of the designer, and/or any guidelines in use and/or any 
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statistical analysis methods that may be employed. Fxirthermore, statistical tools may 
be employed to determine the size of the sample set. 

[0095] With maximum principal stresses determined as a function of strain 
rate and temperature, these stresses can later be used to predict brittle failure in impact 
events. 

IV. USE OF DEFORMATION MODEL AND FAILURE CRITERIA. IN FINITE 
ELEMENT ANALYSES TO PREDICT PART FAILURE: 

[0096] Failure criteria, both ductile and brittle, may be compared to equivalent 
plastic strain and maximum principal stress levels from finite element analyses of the 
part, manually, by looking at result text listings or by looking at contour plots. If these 
levels are above the failure criteria values obtained, then failure may be predicted. 
Preferably, the user defined subroutine of FIG. 15 (which includes the deformation 
model and failure criteria for the material described earlier) is used with a 
commercially available finite element analysis package to automatically compare 
equivalent plastic strain levels and maximum principal stress levels versus ductile and 
brittle failure criteria (which have been predetermined by the method described 
previously). If either failure criteria is exceeded, the user subroutine automatically 
sets the element stiffiiess matrix to zero (wherever the criteria is locally exceeded), 
simulating failure at that location. The load carrying capability of the part will 
decrease as more element stiffiiess matrices go to zero as a resuh of the failure criteria 
being exceeded. 

[0097] A testing and material modeling methodology has been presented to 
model the deformation and failure behavior of painted engineering thermoplastic 
materials. The deformation of the material is characterized by a 5 material parameter 
constitutive model along with a table of total stress minus yield stress versus plastic 
strain. The constants and post yield table are obtained by performing tensile and 
compression tests. The ductile / brittle behavior of the painted material is 
characterized by performing painted disk impact tests. From these tests the failure 
mode that would be expected can be mapped out. Failure criteria are obtained by 
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performing finite element analyses using the 5 constants and post yield table discussed 
earlier in a finite element user material subroutine as shown in FIG. 14. An 
alternative constitutive model for deformation could be used, but the model described 
above is preferred. By correlating stress and strain levels in these finite element 
simulations with experimental failure displacements, failure criteria may be obtained. 
For ductile failures an equivalent plastic strain failure criterion may be used. For 
brittle failures a maximum principal stress failure criterion may be used. Once the 
failure criterion have been established, part performance may be predicted through a 
fmite element analysis using a constitutive model that includes the same deformation 
model used to obtain the failure criterion along with the failure criterion that were 
established. An explicit finite element user material subroutine using the preferred 5 
material parameter and post yield table deformation constitutive model and the ductile 
and brittle failure criteria is shown in FIG. 15. 

[0098] The embodiments described herein account for paint effects upon 
failure behavior, accoimt for biaxial stress states, allow for different failure 
mechanisms and modes and appropriately assign different failure criteria to different 
failure mechanisms. Biaxial painted disk impact tests are used to map out the 
potential failure modes and to generate data for developing failiu-e criteria. Stresses 
and strains at failure are accurately determined by performing finite element 
simulations of the simple test geometries to determine failure criteria (rather than 
attempting to measure a strain to failure value in a tensile test). 

[0099] Note that prior art techniques fail to distinguish between failure modes 
and merely rely on a tensile failure strain to predict failure, thus, not taking into 
account the effect of different stress states. Further, such prior art techniques rely on 
uniaxial stress states in tensile specimens to generate failure criteria. Such prior art 
techniques are deficient because they fail to test biaxial stress states which are usually 
encountered in painted parts. In addition, the effects of the paint upon the failure 
prediction are often not considered. A proven method for accoxmting for the effects of 
paint upon failure behavior in a predictive sense has not been demonstrated in the 
prior art. In contrast, the embodiments described herein account for the effect of a 
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biaxial stress state and account for the effect of the brittle paint system on the ductile 
thermoplastic substrate. Biaxial painted disk impact tests are performed. Failure 
modes are examined as a function of strain rate, temperature, and paint system and 
failure criteria are generated for each condition by correlating coupon test results with 
finite element analyses employing the deformation model to determine stresses and 
strains at failure. In addition to providing a good predictive capability, the approach is 
also practical, using simple tests and commercial finite element packages having a 
user defined material model capability. Adding in simplicity and practical 
application, the paint and substrate materials are treated as a system rather than 
individually. This approach eliminates the added complexity in terms of modeling 
and material testing that would be required if the paint and thermoplastic were 
modeled individually. 

[0100] In addition, the embodiments described herein allow for a more 
accurate determination of ductile failure strains than is currently possible in the simple 
tensile tests used in current techniques. Tests may be performed on coupon 
specimens, and failure loads/displacements correlated with finite element analyses 
using the deformation model to accurately determine true stresses and strains at 
failure. As mentioned, knowing the potential failxire modes of a material, along with 
having accurate failure criteria for the different failure modes, the impact performance 
of the painted part can be more accurately predicted. Further, knowing whether or not 
a failure will occur, the failure mode, and the load or displacement at failure can be 
predetermined. Thus, the number of part testing trials and design iterations required 
to achieve a satisfactory design may be reduced. 

[0101] The description applying the above embodiments is merely illustrative. 
As described above, embodiments in the form of computer-implemented processes 
and apparatuses for practicing those processes may be included. Also included may 
be embodiments in the form of computer program code containing instructions 
embodied in tangible media, such as floppy diskettes, CD-ROMs, hard drives, or any 
other computer-readable storage medium, wherein, when the computer program code 
is loaded into and executed by a computer, the computer becomes an apparatus for 
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practicing the invention. Also included may be embodiments in the form of computer 
program code, for example, whether stored in a storage medium, loaded into and/or 
executed by a computer, or as a data signal transmitted, whether a modulated carrier 
wave or not, over some transmission medium, such as over electrical wiring or 
cabling, through fiber optics, or via electromagnetic radiation, wherein, when the 
computer program code is loaded into and executed by a computer, the computer 
becomes an apparatus for practicing the invention. When implemented on a general- 
purpose microprocessor, the computer program code segments configure the 
microprocessor to create specific logic circuits. 

[0102] While the invention has been described with reference to exemplary 
embodiments, it will be understood by those skilled in the art that various changes 
may be made and equivalents may be substituted for elements thereof without 
departing from the scope of the invention. In addition, many modifications may be 
made to adapt a particular situation or material to the teachings of the invention 
without departing from the essential scope thereof Therefore, it is intended that the 
invention not be limited to the particular embodiments disclosed for carrying out this 
invention, but that the invention will mclude all embodiments fallmg within the scope 
of the appended claims. 
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